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Abstract 



We show that in QCD in the leading twist approximation flavor singlet light 
cone current-current correlators increase with distances rapidly and without 
spatial oscillations. 
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I. INTRODUCTION 



It has been suggested that increase with energy of cross sections of hard processes found 
in QCD is equivalent in coordinate space to the increase with distance of flavor independent 
correlators of currents on light cone ref. [1] . A number of simplifying approximations made 
in ref. [1] like double logarithmic approximation with fixed coupling constant precludes 
unambigous conclusions. In this paper we show that increase of correlators of currents 
with distance is valid in the DGLAP, BFKL and resummation approximations. Space-time 
oscilations of correlators of currents which are present in the parton model [2], disappear in 
QCD for sufficiently large invariant length py. 

Increase with energy of cross sections of hard processes can not continue forever because 
it contradicts to conservation of probability. [3] Application of LO DGLAP approximation 
to the interactions of colorless gluon dipole gg with the virtuality ~ lOGeV^ ^ Ql 
with a nucleon target reveals this conflict with probability conservation in the kinematics of 
X = (Q^/W^ + Q^) < XcriQ^) ^ 10-5 achievable at LHC. [3,5]. (Q^ ^ iQGeV^ was chosen to 
guarantee smallness of running coupling constant.) At the same time the amplitude of the 
scattering of qq dipole in this kinematics is far from S matrix unitarity limit. Onset of new 
pQCD regime at ~ 1 — 2GeV^ for the interaction of qq dipole has been proposed in [6,7] 
and references therein. 

Account of unitarity of S matrix, of energy-momentum conservation law and nonconser- 
vation of bare particles helps to show that limiting behavior of hard processes at sufficiently 
high energies is described by black disc with the radius increasing with energy=BDL regime. 
[4,5]. A lot of new hard phenomena which may appear effective tool to identify BDL regime 
can be observed at LHC . For the review and references see [5]. The aim of this paper is to 
demonstrate that color and flavor blind correlators of currents calculated within conventional 
approximations rapidly increase with relative distance. Such a behavior of correlators is due 
to Lorentz dilatation of life time of colorless dipole produced by current in the target rest 
frame. It resembles turbulence where correlators of velocities are increasing with relative 
distances [8]. Increase with distance of the correlators calculated in LT approximation and 
even faster increase with distances of HT effects suggests similarity between the transition 
to BDL and phase transition. Really it is well known that key property of a system near 
critical point for a second order phase transition is increase of correlators with distance [9]. 
Since our interest is in the dependence of commutator of currents on distance but not its 
absolute value we will not distinguish below between different correlators (different structure 
functions) . 

For dcfinitness we begin our discussion from the DIS lepton scattering off hadron target 
T: e + T ^ e + X. In this case is the virtuality of photon and W is invariant energy 
of 7*T collision. Our interest is in the regime of small x = Q'^/W'^ + which is however 
significantly larger than XcriQ"^) characteristic for BDL. It is well known that total cross 
section of DIS can be calculated as matrix element of correlator of electromagnetic currents 
at light-cone . Thus knowledge of structure functions as a function of x, is sufficient 
to reconstruct correlators of e.m. currents in coordinate space. B.Joffe [2] was the first 
to extract this correlator from experimental data on structure functions of DIS within the 
framework of parton model. In ref. [11] Fourier transform into coordinate space of variety 
of expressions for structure functions has been calculated.In this paper we restrict ourselves 
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by evaluation of commutator of currents at large invariant distances (py) — > oo and small 
but fixed y"^. (Here y is relative space-time difference between currents within a correlator.) 
In this kinematics correlator is dominated by structure function in the kinematics a; -C 1 
and small but fixed Q^. 

The method of moments or Wilson operator expansion for the product of currents is 
convenient to evaluate the behavior of the current-current commutator in the coordinate 
space. 

< N\j^iy)3u{o)\N > = < 7V|O„(0)|7V > +NLT terms 

n 

= P^Pu{py)f{py, y^)/{y''Y + NLT terms 

(1.1) 

where (py) is a kinematical factor. In this paper we shall work with 

G{py, y') = J d^qGix, q') eMm) (1-2) 

where G{x, Q^) is a structure function. 

Another useful quantity is the Fourier transform of gluon density which can be interpreted 
as oc effective coloress "potential" D of dipole-target interaction. 

D{py,y ) = J d q exp(igy) (1.3) 

This "potential" is a Fourier transform of a dipole-target cross-section, and xG is a cor- 
responding gluon structure function. (We use here conventional notations for structure 
functions. The name "potential" is used because within the nonrelativistic approximation 
D has meaning of potential.) D is often used in the eikonal models for the structure functions 
in the regime of small x as "potential" [6,7]) 

We have also studied the closely related problem of coherence length in QCD. We explain 
that coherence length which measures longitudinal length of virtual photon wave function in 
the target rest frame is Ic ~ l/2mNX , far from BDL. However near BDL coherence length 
becomes oc x^~'^ <^ l/2mNX because of a more rapid increase with energy of the interaction 
of smaller dipole with a target, /j, ~ 0.2 as known from theoretical and experimental investi- 
gation of xGp{x, Q"^). This is because of the enhancement of contribution of small transverse 
size configurations in the wave function of photon whose interaction increases with energy 
faster than for average configurations. 

Our results suggest that near BDL the dipole wave function should be modified due to 
the long range dipole-target interactions discussed above. Significant change of the wave 
function occurs in the kinematics where leading twist approximation is violated. 

The paper is organized in the following way. In section 2 we consider the correlators of 
currents in the leading log DGLAP approximation. Then we explain that similar behaviour 
arises in other pQCD approaches also. In section 3 we discuss the coherence length in QCD 
far and near BDL. In section 4 we discuss the deformation of the dipole wave function. Our 
conclusions are in the section 5. 
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II. CURRENT-CURRENT CORRELATORS IN THE LEADING LOG DGLAP 

APPROXIMATION. 



The DGLAP evolution equation in the leading log approximation is 

Q'^.Gix, Q') = (as/{2n)) j\dx' /x'hG{x/x')G{x', Q'). (2.1) 

Here 7g is the kernel in the QCD evolution equation. 

Consider first the case of a freezed coupling constant (considered in ref. [1]). In this case 
the structure function is given by 

G(x, Q') = I j^,x'W(j, as log(gVgg)) (2.2) 
Jc [Zm ) 

The contour of integration over j (we denote it C) runs along a straight line parallel to the 
imaginary axis to the right of all singularities of the integral. We use the notation = —q^ 
if < and = if > 0. The Bjorken scaling variable is defined in a usual way: 

x^-qyi2pq). (2.3) 

For the anomalous dimension we have: 

V±{j) = OMVfU) + Vg{j) ± ^/{Vf{j) - VG{j)r + 24$G(^-)$g(j-), (2.4) 

where 

VfU) = -C2(4V^(j + 1) + 47^ - 3 - 2/{j{j + 1)), (2.5) 



Vg{j) = -4iVe(V'(j + 1) + 7e) + lliVc/3 - 2 + 8N,{f + j + l)/{j{f - + 2)), (2.6) 



$^(j) = 2C2{f+j + 2)/{j{f - 1)), (2.7) 

*S(j) = (/+J+2)/(j(j + l)(j + 2)). (2.8) 
The function is the usual logarithmic derivative of Gamma function 

and 

7g(j) = V+U) (2.9) 

Let us begin calculations within the simplifying assumption of frozen coupling constant 
and then to analyse realistic case. It is easy to carry out the integration over d'^q in the 
Fourier transform. We use the formulae [1]: 
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{qyF{x) ^ 7rr(/3 + 2)M^ F{x)x'' cos{x{py) + (5'K/2)dx (2.10) 

if the latter integral converges. Here w+ = 9{u)u. For the case of the power hke behavior 
of F{x) this integration can be performed analytically, and in the limit of large invariant 
length {py) ^ 1 one obtains power like asymptotics: 

{qyx^ ^ nr{P + 2)r{a + P + 1) cos((2a + P)n/2)/{{pyr){yy+') (2.11) 

In our case P — 7(j), a — —j, so 

G{y\mj) = (V2) / dj{mj)^'-^^^^^n^G{j) + 2)r(l - j + ^aij)) sin((-j + 2^G{j))^m 

(2.12) 

We now use the saddle point method to integrate over j. 

log(l/oVy^) d^cU) 
log(py) dj 

The saddle point can be found numerically. The derivative of the DGLAP anomalous di- 
mension can be approximated with good accuracy by 

i{n) ~ -l2/{n-lf 

for 1 < n < 2.5 and 

y ~ -5/(n- 1) 

for n > 3. This means that the double log expression works quite good numerically for 

\og{py)>as\og{yl/y^) (2.14) 

Using the latter numerical approximations for anomalous dimension, one obtains the relevant 
saddle point: 



3^1 + Vl2aslog(?/oV?/')/log(p?/) (2.15) 

i.e. we have the double log expression. However now we are permitted to take into account 
also the terms that contain only one big logarithm. For \og{py) > log(l/|/^), one has 

^^Py^ = l!^^ W Wi^<^sNMlog{py)log{Qly-)] 

X cos{{,J{asNjn){^log{Qly^)/log{py) + 2^log{py)/ log{Qly^))n/2} 
X r{^(aM^){^^og(Qly^)/\og(py) + ^Jlog(py) / \og(Qly^))} 



X r{2 + ^J{asNJn){log{py)/log{Qly^)} 

(2.16) 
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Consider now the realistic case of the running couphng constant. In this case [15,17,16] 
the gluon structure function is given by 

Gix,Q')= I ^x-W{j,0 (2.17) 
Jc [Zm ) 



Here 



where 



and 



e = ^logK(gg)/a,(g2)), (2.19) 



6 = 11- -Np = 9 , Nf^3 (2.20) 

3 

The calculations one needs to carry through for the realistic case of the running structure 
constant are completely similar to the calculations with the freezed coupling constant. The 
only change in the calculation procedure is the substitution in all calculations of the terms 
like log((5^/(5o) by log(yQ/y^). Consequently, we have to substitute ccg log(yQ/y^) by 

ay') = llog{as{W)Ml/yl)). (2.21) 
The saddle point will be in the point 



and we obtain 



J = 1 + A -j — 7 — r (2.22) 

\ log(p7/) 



X cos((Ve(y^)/log(p|/)) X r(((^e(y^)/log(p|/))) 



(2.23) 



Similar calculations can be carried for an opposite case \og{py) < ag logiyl / y"^) , ^^"^ one can 
easily see that the correlator increases with the increase of the invariant length in this case 
also. 

Similarly it is easy to derive the increase with distance of the effective dipole-target 

potential. 

One of the important features of the correlators in the parton model is their space- 
time oscillations [2]. We will show that oscillations at large distances in the correlator and 
"potential" for dipole -target interactions disappear in pQCD. 
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For fixed and frozen coupling constant the correlators and the effective dipole-target 
potentials derived above oscillate, as shows the existence of the cos term. The period of 
oscillations can be determined from the condition that the argument of the cosinus is m: 
where n is the integer number. 

7(jo)+Jo = ri (2.24) 

where jo is the saddle point. We immediately obtain the condition 

^a,7V,/7r(^log(l/y2)/log(pt/) + ^log(py)/ log(l/y2)) = n (2.25) 

In our approximation the second term is bigger than the first and we obtain that the oscil- 
lation period is 

py - (l/^2^)nVKiVe/7r) ^2.26) 

Thus oscillations exist in the case of frozen coupling constant for large py. This analysis can 
be easily upgrated to analyze the case of the running coupling constant. In this case the 
relevant condition is 

e(y^)/log(py)=n^ (2.27) 

Hence 

py ~ exp(e(y')/n2) (2.28) 

Thus taking into account of running coupling constant changes dramatically the oscillation 
period: instead of increasing with n it now decreases. Therefore oscillations are present at 
small (py) which are beyond of region of applicability of our consideration. 

The reason for such dramatic infiuence of asymptotic freedom is evident in the procedure 
used to calculate the integral over = —Q^. This integral tis the product of rapidly 
oscillating cos(g^(t — r)/xB) and -F(g^) [2] 



G{px,x^) = - 2^(^py^ J J dxBG{xB,Q'^){q^/x%)s\n{q^/{w?XB)x'^/ipy) + {py)xB) 

(2.29) 



In the case of frozen coupling constant the function G behaves as the power of Q^, and 
rapidly changes within the period of oscillations of cosinus. Thus exact calculation of the 
integral produces terms proportional to cos(y(j)7r/2, or corresponding sinus because of 
F{Q^) ~ Q'^^^^\ This gives rise to the above result . However in the case of the running 
couphng constant the function F{Q'^) depends on Q"^ more slowly , due to the presence of 
the additional logarithm. Indeed, the characteristic period of oscillations is ~ x{py)/y^. 
However the characteristic x is connected to characteristic (py) (coherence length) as py ~ 
l/(x)^ where /3 is a number close but smaller than one. Then the period oscillations is 
~ {pyY~^^^ /y"^- Thus the period of oscillations decreases with increasing [py). Evidently , 
the function Oy'^) on this scales can be considered constant, and one can safely carry out 
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$,{q^) from the integral and substitute it by ^{y^)- The oscillations die down for running 
coupling constant and large {py). The behavior for the case of \og{j)y) < Q;slog(l/y^) 
is qualitatively similar. Thus the oscillations we found within the DGLAP approach are 
qualitatively different from the oscillations in the parton model found in ref . [2] . 

Similar analysis can be carried through for the BFKL case [18] by substituting Qo — 
l/(t — z) in the Fourier transform as it is legitimate in the leading ag ln(xo/x) approximation. 

Formally, we use for the BFKL case the contour integral representation Let us remind 
that in the BFKL case the structure function G can be represented as 

G{x, Q^) = J (dM/27r) exp(Mt + a,x{M)u) (2.30) 

where t = \og{Q'^ /Ql),u = log(l/a;). In the BFKL case the coupling constant is freezed. 
The integration contour is defined in the same way as in the DGLAP case. For the Regge 
limit the integral can be taken using saddle point method. The result is the energy indepen- 
dent saddle point Mq = 1/2. It is easy to take the Fourier transform and obtain according 
to the formulae above, 

G{py, y^) = ^ exp(a,x(l/2) log(py)) (2.31) 
Here the BFKL anomalous dimension x is given by 

^^M) = —(2^(1) - ^(M) - - M)) (2.32) 

TT 

with 

x(l/2) = -log(2) (2.33) 

TT 

We obtain that the correlators are quickly increasing with the distance and the oscillations 
are absent. 

Within the resummation models of the ABF [12] and Ciafaloni et al [13] it was argued 
that the behavior of the structure functions in these models for the wide range of x is the 
same as in leading+plus nonleading order DGLAP. Numerically this behavior is qualitatively 
quite similar to the leading order DGLAP, except that the splitting functions is 15-30 % 
lower than the leading order DGLAP, in the wide range of Bjorken x. At quite low x (say, 
X ~ 10"^ for Q"^ ~ 4.5 GeV), the DGLAP behavior of the kernel switches to BFKL type 

rise. (Note [14] that the switch point depends on very weakly, x^ ~ 1/ ^Jas{Q^)). This 
rise, although its analytic form is unknown leads to the rapid rise of the correlator with 
distance without oscillations 

G ~ ipyr^-' (2.34) 

Here ap is the renormalized perturbative Pomeron slope of the resummation models, which 
is quite close (effectively ap ~ 1.25). Note that the oscillations are absent starting from 
sufficiently large py already in the DGLAP curve. 



7 



III. COHERENCE LENGTH IN QCD 



Let us now evaluate coherence length Ic for the total cross section of DIS evaluated 
within DGLAP aproximation. We use here the same definition as B.Joffe length within 
parton model [2,19] and perform the inverse Fourier transform: 

/■oo roc 

xG(x, Q^) = 2u / dy^ / duF(y^, u) cosivy"^ / i2u) - xu), (3.1) 
Jo Jo 

(the latter equation defining function F [19]). Let us carry on firstly the integration over 
y^. The relevant scale for the integration over is l/Q^- Thus within the LO DGLAP 
approximation: one may substitute y^ ^ so that the integral takes the form: 

roo 

xG(x,Q^)^ 4:uduF(l/Q^,u)sm(l/2xu-xu) (3.2) 
Jo 

The coherence length is defined for given Q^,x as the typical; u in the integral. Such u is 
either determined by the maximum of the function F, or by the radius of oscillations in the 
sinus multiplier. (For the another definition of Ic and therefore different dependence on x 
see Y.Kovchegov and M.Strikman [20].) 

Ic within a parton model was determined by B.Ioffe who assumed that 

F{u) const. (3.3) 

for u oo. If so essential u in the integral are completely determined by the oscillating 
multiplier , i.e. 

ui{x,Q^)^l/x (3.4) 

The same behaviour is expected within DGLAP, LO+NLO BFKL, resummation approxi- 
mations. This is because in all approaches Ic is effectively fife time of dipole with mass M 
where = k^/z{l — z). Here kt is transverse momentum of quark within dipole and z is 
the fraction of dipole momentum carried by quark. Thus Ic ~ Qo/Q^ + and ~ M^. 
Therefore Ic ~ l/2mjsix . 

Near unitarity limit fast increase with energy of the coherence length should somewhat 
slow down. This is because the fast increase of amplitude with energy leads to the related 
increase of probability of configurations of constituents with large kt. Let us consider for 
definitness the case of longitudinally polarized photon where maximal important kt can be 
evaluated from pQCD formulae as: Q'^/k'^l/x^^^*^ ~ 1 So near unitarity limit: Ic ~ Qo/M^ cx 
<^ l/2mNX. 

IV. ON THE DISTORSION OF DIPOLE WAVE FUNCTION NEAR UNITARITY 

LIMIT. 

Another consequence of the increase of the correlators with the invariant length py is the 
deformation of the wave function of the dipole near BDL. 
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We found that fast increase with energy of LT contribution transforms in coordinate 
space into fast increase with distance of correlators of currents . The contribution of higher 
twist effects is increasing with x decrease even faster :cx l/x^^^^^ ^ where n=l for LT n=2 for 
first HT etc. Performing the same Fourier transform of this formulae into coordinate space 
as for LT term will produce even faster increase with distance of the contribution of HT 
effects as oc [pyY+'^i^iQ ). Thus dipolc-targct interaction near unitarity limit is characterized 
by important role of long range interactions of dipole with target (which however are sig- 
nificantly smaller than B.Joffe length). This long-range target-dipole interaction effectively 
leads to the interaction between the components of the dipole. Indeed, at moderately small 
X the contribution of dipole configurations with z ~ 1/2 is multiplied in the cross-section 
by the small factor -square of the dipole moment d"^. Near BDL this suppression disappears 
and the dominant configurations are z ~ 1/2. Thus effective wave function of the dipole 
changes significantly near the unitarity limit. 

V. CONCLUSION 

The main aim of the present paper was to demonstrate that the correlators of currents 
fastly increase with distance in pQCD . 

We have also studied the issue of the space-time oscillations of the correlator of currents 
and found no oscillations in QCD contrary to the parton model where the correlator of cur- 
rents is ~ cos{py). Within DGLAP approximation oscillations are decreasing with distance 
because of decrease of coupling constant with virtuality. 

We estimate coherent length near BDL as Ic oc a;^"'^ . 

We argue existence of the deformation of the dipole wave function near BDL resulting 
from the long range interaction between dipole and target [1]. 

The most suggesting consequence of our results is that in condensed matter physics 
increase of correlations with distance usually leads to the phase transition. [8] 

One of us (LF) is indebted to M.Strikman for the discussion of increase of correlators 
with distance which led to the understanding of the impact of increase of gluon distribution 
on the coherence length. 
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